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• Premise of the study: The genus Ficus (fig trees) comprises ca. 750 species of trees, vines, and stranglers found in tropical 
forests throughout the world. Fig trees are keystone species in many tropical forests, and their relationship with host-specific 
wasp pollinators has received much attention, although many questions remain unresolved regarding the levels of host specific- 
ity, cospeciation, and the role of hybridization in fig and wasp speciation. We developed exon-primed intron-crossing (EPIC) 
markers to obtain phylogenetic resolution needed to address these questions. 

• Methods and Results: Expressed sequence tags (ESTs) from F. elastica were compared to Arabidopsis and Populus genomes 
to locate introns and to design primers in flanking exons. Primer pairs for 80 EPIC markers were tested in samples from diver- 
gent clades within Ficus and the outgroup Poulsenia (Moraceae). 

• Conclusions: Thirty-one EPIC markers were successfully sequenced across Ficus, and 29 of the markers also amplified in 
Poulsenia, indicating broad transferability within Moraceae. All of the EPIC markers were polymorphic and showed levels of 
polymorphism similar to that of the widely used internal transcribed spacer (ITS). 
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Ficus L. (Moraceae) is a pantropical genus comprised of 
ca. 750 species of trees, epiphytes, shrubs, vines, and stranglers 
found primarily in humid tropical forests. As a year-round source 
of calcium-rich fig fruits, Ficus trees are often described as 
keystone species. However, Ficus may be best known for their 
pollination mutualism with small (1-2 mm), short-lived (1-2 d) 
"fig wasps" in the family Agaonidae (Weiblen, 2002; Herre et al., 
2008). Female fig wasps pollinate flowers and oviposit within the 
enclosed inflorescence (syconium or "fig"), in which the larvae 
develop before emerging to pollinate and oviposit in the syconia 
of asynchronously flowering conspecific trees. For sustained 
reproduction of the figs and the wasps, the wasps must exhibit a 
high degree of host- specificity, and the host population must pro- 
vide access to flowers (i.e., figs) throughout the year. 

Although the fig-wasp pollination mutualism is one of the 
tightest known in terms of host-pollinator specificity, there are 
many exceptions to the one pollinator species/one host species 
rule. In some cases, two or more wasp species pollinate the 
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same host species in different parts of its geographic range, and 
multiple wasp species have been found in a single host tree 
(Herre et al., 2008). Furthermore, in Central America and in 
South Africa some wasp species have been shown to use more 
than one fig species in the local fig community (reviewed in 
Herre et al., 2008). The nonspecificity of some pollinators, in 
addition to some genetic studies (e.g., Machado et al., 2005), 
suggests that hybridization is possible. 

Most phylogenetic studies of Ficus have used chloroplast 
DNA and/or one or two commonly used nuclear DNA markers 
(e.g., internal transcribed spacer [ITS]) (e.g., R0nsted et al., 
2005). These markers are insufficient in number for studies of 
introgression, and they do not resolve phylogenies of closely 
related species or phylogeographic structure in widespread spe- 
cies (C. Dick, unpublished). To address the deficiency in nuclear 
genomic markers for Ficus, we have developed a set of exon- 
primed intron-crossing (EPIC) markers by comparing an ex- 
pressed sequence tag (EST)-library for F. elastica Roxb. ex 
Hornem. with the annotated genomes of Populus trichocarpa 
Torr. & A. Gray (Salicaceae) and Arabidopsis thaliana (L.) Heynh. 
(Brassicaceae) using a bioinformatics pipeline developed by Li 
et al. (2010). 

METHODS AND RESULTS 

Selection oftaxa — Neotropical Ficus contains two distinct and phylogenetically 
distant subgenera, which represent two important neotropical life forms: the 
free-standing fig trees (subg. Pharmacosycea (Miq.) Miq. sect. Pharmacosycea) 
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and the strangler figs (subg. Urostigma (Gasp.) Miq. sect. Americana Miq.). 
Sect. Pharmacocysea is sister to all the other fig subgenera, and therefore our 
sect. Americana and sect. Pharmacocysea samples share a most recent common 
ancestor that is the base of the entire Ficus crown clade, which, based on fossil 
records, dates back to at least 60 million years before present (R0nsted et al., 
2005). All primers were tested on F obtusifolia Kunth (sect. Pharmacocysea) 
and F maxima Mill. (sect. Americana), which were collected from the Barro 
Colorado National Monument (BCNM) in central Panama. The subset of 
primers that amplified in both Ficus species were also tested on Poulsenia ar- 
mata (Miq.) Standi., which is a monotypic genus in the fig family Moraceae 
(Datwyler and Weiblen, 2004). Botanical vouchers (Dick and Gomez 234, 
F. obtusifolia; Dick and Gomez 240, F. maxima; and Dick and Gomez 180, 
P. armata) were deposited at the herbaria of the University of Panama (PMA) 
and University of Michigan, Ann Arbor (MICH). Genomic DNA was extracted 
with the cetyltrimethylammonium bromide (CTAB) method of Doyle and 
Doyle (1987). 

Bioinformatics pipeline — Researchers from the United States Department 
of Agriculture (USDA) previously developed an EST library of F. elastica to 
characterize the genetic basis of rubber biosynthesis (McMahan and Whalen, 
personal communication). We compared 9289 unique F. elastica ESTs from 
the National Center for Biotechnology Information (NCBI) database with the 
annotated genomes of A. thaliana (Brassicaeae) and P. trichocarpa (Salicaceae) 
using the informatics pipeline developed by Li et al. (2010). Briefly, we (1) re- 
trieved coding sequences (CDS) that were longer than 100 bp from the anno- 
tated genomes of A. thaliana and P. trichocarpa. (2) We compared those CDS 
with the genome of the same species to identify "single-copy" CDS. (3) The 
candidate single-copy CDS thus identified were subsequently compared to the 
EST library of F. elastica to find markers that were conserved (identity >80%) 
among all three species. (4) After locating the single-copy conserved CDS, we 
screened for CDS flanking small introns, which were smaller than 1000 bp in 
the compared genomes, to facilitate the subsequent PCR and sequencing steps. 
Primers based on the F. elastica exons were initially designed by eye and 
subsequently checked with the Primer3 web interface program (Rozen and 
Skaletsky, 2000). 

Primer assays — PCR was performed in a final volume of 20 jiL containing 
10 mM Tris-HCl (pH 8.4), 50 mM (NH 4 ) 2 S0 4 , 1.5 mM MgCl 2 , 0.2 mM dNTPs, 
0.1 jiM each primer, 2 ng of genomic DNA, and 0.5 units of Taq polymerase 
(BioTherm, Gaithersburg, Maryland, USA). The amplification profiles included 
an initial denaturing at 94°C for 5 min; followed by 35 cycles of 50 s at 94°C, 
50 s at 54°C, and 1 min at 72°C; and a final extension step of 10 min at 72°C. 
PCR products were ligated into the pMD 18-T plasmid vector (Promega Corpo- 
ration, Madison, Wisconsin, USA) and transformed into E. coli strain (DH5a, 
Promega Corporation). Insert-positive plasmids were isolated using the 
E.Z.N.A. Plasmid Mini Kit I (Omega Bio-Tek, Norcross, Georgia, USA) and 
amplified using M13 primers. Forward and reverse strands of each amplicon 
were sequenced on an ABI 3730xL DNA sequencer (Applied Biosystems, 
Carlsbad, California, USA) at the University of Michigan Sequencing Core 
Facility. All Ficus insert sequences have been deposited in GenBank (accession 
numbers JQ341915-JQ341980; also see Table 1). For comparisons with ITS, 
we also obtained ITS sequences for F. obtusifolia, F. maxima, and P. armata 
(GenBank accessions JX1371 13-JX137114) using standard methods. 

Data analyses — DNA chromatograms were edited using the SEQUENCHER 
program (Gene Codes Corporation, Ann Arbor, Michigan, USA). DNA sequences 
were initially aligned using ClustalX version 1.81 (Thompson et al., 1997) with 
default settings, and subsequently aligned manually using Se-Al (Rambaut, 
1996). We determined number of polymorphic sites, nucleotide diversity (ri), 
and GC content using MEGA 5 software (Kumar et al., 2008). 

Results — We identified 200 ESTs that satisfied our criterion of 80% exon 
identity with the published genomes. Based on intron length, we selected a 
subset of 80 ESTs for further marker development, of which 31 amplified success- 
fully in Ficus species from both subgenera, 16 amplified in one species only, 
and 33 did not amplify in either species. The 31 cross-amplifying primer pairs 
were further tested in P. armata, of which 29 amplified successfully (Table 1). 
The number of polymorphic sites in F. obtusifolia and F. maxima comparisons 
ranged from 12 to 71 (mean = 32), whereas nucleotide diversity ranged from 
0.02655 to 0.07305 (mean = 0.0470) (Table 1). In comparison, there were 45 
variable sites in ITS between F. obtusifolia and F. maxima, falling within the 
range of the EPIC marker variation. 
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CONCLUSIONS 

The 31 EPIC markers that amplified between the two Ficus 
subgenera indicate that these markers might be useful across the 
full phylogenetic breadth of the >60 Ma genus and its >750 spe- 
cies. The markers that transfer to Poulsenia indicate an even 
broader phylogenetic utility within the Moraceae (ca. 40 genera 
and 1000 species), which probably originated in the Cretaceous. 
These markers should therefore be extremely useful for phyloge- 
netic analysis at the family level and potentially beyond. The 
markers show a level of intron divergence that is of a similar 
magnitude as ITS, which is one of the most informative and 
broadly used markers in plant molecular systematics. These EPIC 
loci should be useful for analyzing recent divergences in which 
incomplete lineage sorting and/or introgression may be factors, 
including recent speciation, hybridization, and comparative phy- 
logeography. In combination with EPIC markers developed for 
chalcid wasps (Lohse et al., 2011), it should now be possible to 
jointly analyze wasp and host plant phylogenies to study revo- 
lution at both population and phylogenetic scales. 
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